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I. 


INTRODUCTION 


This is a study of the feasibility of obtaining optical images of a 
cometary nucleus via a flyby of Comet Encke. This study is based on a 
physical model of the dust cloud surrounding the nucleus. Development of the 
model is based on available physical data and theoretical knowledge of com- 
etary physics. Using this model, calculations are made of the absolute 
surface brightness of the dust in the line of sight of the on-board camera 
and the relative surface brightness of the dust compared to the nucleus. 

The brightness is calculated as a function of heliocentric distance and for . 
different phase angles (sun- comet -spacecraft angle). 

The .study of the feasibility of obtaining optical images of the nucleus 
of Comet Encke via a flyby was divided into two parts. First a physical 
model of the 'dust cloud surrounding the nucleus was derived and then our 
Mie scattering code was employed to calculate the absolute brightness 
of the dust in the line of sight of the on-board camera. The surface 
brightness of the comet was compared to this background brightness along 
with contributions' from debris and ice particles. Conclusions were then 
drawn .in light of these results regarding the optimum imaging system for 
a flyby mission. 

H. BACKGROUND OF FLYBY MISSIONS 

At a symposium on the exploration of space held in Washington, D.C. 
in April 1959 > Whipple (1959) pointed to the possibility of sending a space 
probe through the neighborhood of a comet. Two years later. Swings (1962) 
gave an extended survey on the scientific objectives and feasibility of 
such a mission at a symposium held in Pasadena, Calif, in August 1961. 
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He summarized the existing knowledge about comets and pointed to those 
X>roblems in the physics of the comets for which investigations by means of 
space probes would be of special value. 

There are quite a few problems which would be handled best with instru- 
ments installed in a cometary probe (Lust, 19 69 ) : 

(a) The structure and composition of the nucleus , its surface tem- 
perature and its color. 

(b) Chemical processes related to the formation of the observed radi- 
cals which take place in the vicinity of the nucleus where the density is 
high, and the ejection velocity of the different particles. 

(c) The composition of the cometary atmosphere, the mechanisms of 
dissociation and ionization, the size and nature of the dust grains. 

(d) The density gradient and partial densities of the different con- 
stituents in the coma and in the tail. 

(c) The interaction with the solar wind and with interplanetary mag- 
netic fields, the structure of shock fronts and related problems of plasma 
physics and magneto-hydrodynamics - 


Among the questions which have to be solved in the preparation of a 
cometary probe in which comet should be selected for a first space mission. 
Much work has been put into this problem by different groups in the U. S.A. 
and in Europe, and all these groups came to realize that the difficulties 
of such a mission are far greater than had first been anticipated. It is 
evident that a periodic comet which has been observed for many apparitions 
and whose orbit has been calculated with some accuracy is an easier target 
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•than a new comet which- suddenly appears without being predicted and whose 
path hah to be calculated from a few observations taken in small time inter- 
• vale. Bui, there: axe other reasons why a periodic comet can be reached more 
•easily by a space vehicle than a new one. The orbits of almost all periodic 
comets have relatively small inclination to the ecliptic plane, and a 
motion in the same' sense as Jupiter and all other planets (direct orbit). 

The .inclinations of the orbital planes of the new comets are, on the other 
hand, randomly distributed between 0 and 1B0°, so that on the average these 
comets spend only a small fraction of their time close to the ecliptic plane * 
while they axe passing- through the nodes. This causes a severe restriction 
for the launch window, while retrograde motions (inclination between 90 
and l80°) cause very large approach velocities of the spacecraft with 
respect to the comet. Expensive terminal guidance would therefore be 
necessary to force the trajectory of the spacecraft into or near to the 
orbital plane of the new comet aid to diminish the relative velocity to a 
slow fly-by or to achieve a rendezvous. 

Because of these difficulties, the new comets have been excluded as 
a first aim for a cometary probe by all groups who have investigated the 
feasibility of such a mission in spite of their high scientific interest. 


Different groups in the U.S.A. as well as in Europe have carried out 
detailed feasibility studies during the last years. In 1961, a space missidn 
to a comet ms considered for the first time in NASA contracted studies of 
the Scientific and Technical Laboratories Inc. (Corben, 1962). 
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There exist about 50 periodic comets which have been observed during 
inoro than one perihelion passage, and an additional 30 comets with perxods 
less than 20 years observed during one apparition only. About 5-7 ° 1 these 
cornets appear on the average per year. A great amount of work and extensive 
calculations have been carried out by the different study groups to select 
among these comets the objects being fitted for a first mission. Several 
selection criteria had to be applied with respect to the following points: 

(a) The position of the comet at the time of intercept must be known 
with very high accuracy, and it became more and more clear that this condi- 
tion is a major constraint for a cometary mission. It seems desirable that 
the probe should approach within 1000 km. If costly 
midcourse corrections are to be avoided, the position of the comet must 
be predictable at the time of launch, that is several months before 
intercept, with this accuracy. 


To secure the necessary accuracy in the calculation of the orbital 
elements, it is necessary to start with a well-known orbit calculated 
from previous apparitions and to apply corrections by means of new observa- 
tions made after recovery. Perturbations by planets, especially by Jupiter, 
can be the reason for considerable changes of the orbit. Also secular 
perturbations caused by a mass loss of the comet and non-gravltational 
forces play a role in this respect. Especially the location of the comet 
in its orbit, that is the daily motion, is affected by such perturbations. 
Since small errors in the daily motion will add up and result in large errors 
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in the time of perihelion passage, the accuracy of this parameter is very 
important (Porter etal., 1965 ). It seemed, therefore, necessary to exclude 
all comets which had "been observed during one apparition only, and also 
some of the remaining comets whose predicted orbits were not well enough 
established. 

(b) In order to guarantee observations of the comet’s, position 
necessary for orbit corrections, the comet should be .recovered at least 
2 months before launch, and it should be seen against a dark sky for several 
hours per night. This imposes the condition that . the comet must be 
at least of 20th magnitude -- the limit for recovery of an object 
whose position is known -- several months before perihelion passage. 

The average brightness of periodic comets is so .small that not many of 
them fulfill this condition. , 

The question of the accuracy of cometary orbits has been investigated 
by the groups in the U.S.A. as well as the European ESRO group. The conclu- 
sion is that if every effort is made to recover a comet very early, that is 
about 20th magnitude, it should be possible to correct the predicted orbits 
to the accuracy necessary for a space mission for a number of comets. 

(c) The comet should be bright enough at intercept to get photome- 
tric data and good quality spectra from ground-based observatories to 
supplement the spacecraft data. This means that the comet should be of 12th 
magnitude or brighter at the time of perihelion passage, and it should be 
seen against a dark sky for some hours. Since a great percentage of the 
periodic comets never become so bright, and since some of the brighter ones 
are in unfavorable positions on the sky near perihelion ( close to the sun) , 
this is a further restriction of the number of feasible comets. 
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(d) The ideal launch velocity should not exceed ^ 1 6 km/sec. This 
limit is chosen rather arbitrarily, it represents for instance the velocity 
required for a two year flight to Jupiter. It turned out, however, that 
fills requirement does not impose a severe restriction, because almost every 
comet that was brighter than 12th magnitude at intercept fulfilled this 
condition. 


(e) Of two otherwise equal missions, that which leads to a smaller 
relative approach velocity between spacecraft and comet and therefore 
allows for a longer stay of the spacecraft in the vicinity of the comet 
should have priority. 

(f) All the points mentioned before refer more or less to the 
technical side of the mission. When choosing a comet, one should of 
course ask which mission promises the best scientific results. Comets are 
of very different activity according to their "age", that means according 
to the number of perihelion passages they have already made, because with 
every approach to the sun they lose a considerable amount of volatile 
material, and slowly get exhausted and inactive. 

Summarizing the different points, a feasible comet should fulfill 
the following conditions: 

(1) A reliable predicted orbit should be available, 

(2) The comet should be of at least 20th magnitude 2 months before 
launch, and at that time be well observable for at least 2 hours per night. 

(3) It should be of at least 12th magnitude at intercept, and dis- 
tinctly visible from the earth. 
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(4) 'The launch velocity should not exceed ~ 16 km/sec. 

(5) The relative approach velocity between comet and spacecraft 
should be low, 

(6) The comet should be interesting from a scientific point of view. 

A mission to Comet Encke in 1980 meets these requirements and should be 
undertaken. 
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III. PHYSICAL PROPERTIES 

With an estimate for the size of the nucleus of Comet Encke , and 
using the icy- conglomerate model, we estimate the distribution of non-volatile 
dust particles near the nucleus. 


A. Nucleus 

Although Whipple's (1950) model for comet nuclei is not universally 
accepted (e.g., Lyttleton, 1972) > it has been used to successfully explain 
a number of comet ary phenomena. Not the least of these is the non-gravitational 
force required to rex>roduce the orbital motion of Comet Encke and several 
other comets (Maraden and Sekanina, 197 + and earlier papers in this series). 

The non-gravitational acceleration derived by Marsden and Sekanina (197 1 *) 
for Comet Encke shows a regular decrease which is consistent with the existence 
of a porous, rocky core, within which "dirty ice" is embedded (Sekanina, 1969b). 
Since the gravitational acceleration is proportional to the fractional rate 
of change of mass (as well as the degree of anisotropy of the ejection), it 
is necessary to estimate the mass of the nucleus in order to obtain the 


rote of mass loss itself. The mass (or mean density) and radius are also 
required in order to calculate the terminal velocities of ejected dust parti- 
cles,, since expanding vapor from the nucleus must lift the particles against 
gravitational attraction (Delsemme and Miller, 1971). Marsden and Sekanina 
(1971) estimate that Cornet Enke loses 0.03 percent of its mass during each 
orbital revolution. A spherical nucleus of radius = 1.7 ^ and mean den- 
sity 1 g cm' 3 has a mass ~ 2 x 10 l6 g, implying AM ~ 6 x 10 12 g per revolution; 
this mass loss is comparable to the estimates of Sekanina (1969a), which are 
~ lO^g. If the geometric albedo of the nucleus is = 0*1, then 
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P 2 ^ which is to be compared with observed values which range 

n n 

from 0.24 (Roemer, 1966 ) to 7 > 0.82 km 2 , which we derive from Roemer' s (1972) 
observation of the apparent magnitude of Comet Encke at aphelion(m ~ 20.5)- 
These lower 1 limits result from our assumption that the comet was 
observed at opposition, thereby giving minimum earth-comet distance 
and phase- -angles These results are consistent with the recent treatment 

of Delsemme . and Rud (1973). 

B. Non-Volatile Particles 

1. glze Distribution 

Lacking a determination of particle ejection for Comet Encke. itself, 
we use the results of Finson and Probstein ( 1968 ) for Comet Arend-Roland 
and those of Sekanina and Miller (1973) for Comet Bennet. In. particular, 
we use the distribution of particle sizes deduced by Finson and Probstein 
( 1968 ), which they state is well -determined over the range 4 < pd s 40 pm g cm 
where p is the bulk density of a particle of diameter d. We use the Sekanina 
?ind Miller (1973) distribution for the smallest particles (down to pd - 0*9M^ 
g era" 3 ). Specifically, denoting by n(a) da the number of particles with 

radii in the range a to a+da, 


K 3 _( a/ a Q ) 5 (a-l)/a Q 

■ 1 <: a £ 6.7 pm 

Kg(a/a Q ) 3 

6.7 < a ^ 14.4 pm 

K^(a/a o )" 4 

14.4 <a < 44.4 pin 

K 4 (a/a o )' 5 

a > 44.4 pm 
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The distribution for larger particles fits smoothly onto the distribution 

_R 

t,y inf e r plan e tary particles as given by Whipple (1968): n~ a . These 
are shown in Figure 1. The constant a - 1 pm, and the K' s are determined 
by the requirement that 

03 1 

(4itp/3) [ a n(a)da = 1 g. 
a 

o 

The rate of ejection of particles of radius a is then given by 
P(a) = 11 (a) j 

where dM/dt is the total mass-loss from the nucleus in the form of dust. 

We have 

Ik 

K ± = 8.916 x 10 

K 2 = 1.137 * 10 14 
K 3 = 1.642 x 10 15 
K 4 - 7.298 x 10 l6 . 

-3 

The particle mass density is assumed to be p = 0.45 g cm , 

which is approximately the value adopted by Whipple (1967, 1968) for meteroids 

of mass m ^ lO - ^ g. Super-Schmidt data for the Taurid meteor shower give 

3 

a mean value of about 0.28 g cm"“ (Verniani, 1969). A somewhat higher den- 
sity may be appropriate for the smaller particles: Whipple (1968) uses 

p'-lg cm” J for m ^ 10 g (a < 175 pm), and radio data for sporadic and 

-3 

shower meteors give p ^ 0.8 g cm (Verniani, 1973)° The particle distribu- 

-3 

tion would not be greatly changed if we were to use p 1 g cm , except 
that n(a) would peak at a ^ 0.6 pm and extend to about 0.55 pm. These 
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LOG a (CM) 


Figure 1 Particle size distribution. 
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small particle;:, would make only a email contribution to the total scattered 
light. 

It is of interest to compare our distribution to one derived by Taylor 
et al. (1973) in a similar manner. This is done in Fig. 2. The disparities 
arise from the different densities and size cutoffs used. 

2. Production Rate 

As discussed above, non-gravitational acceleration of Comet Encke 
implies a loss of - 10 l3 g per orbital revolution. It is not known how 
much of this is in the form of dust, but Whipple ( 1967 ) suggests that 
the abundance of Taurid meteors is evidence that perhaps an old comet such 
tU! Encke releases a larger fraction of material in dust than does a "new" 
comet (see also Delsemme 1 s 1973 review). Whipple ( 1967 ) goes on to esti- 
mate that, averaged over one period, Comet Enke contributes as much as 
3 5 x 10 6 g sec”^" of meteor! tic material to the interplanetary cloud. If 
true, this would imply a much smaller degree of anisotropy of ejection 
than has been determined for this comet (Marsden and Sekanina, 1971; P* 

The distribution of particle sizes might differ from that for "new" comets, 
having fewer small particles. 


We conservatively assume that the tot; el mass lost per orbit is 
of which 10 percent is dusf. This gives a simple average mass-loss 
rate of dM/dt =*• 9.6 x 10 3 g sec ^ as dust. As an approximation to the 

evaporation rate of water snow (Delsemme and Miller, 1971; Marsden, Sekanina 
and Yeomans, 1973), the mass-loss varies as r 2 for r ^ r £ amd is zero at 
greater distances. The mass-loss rate at any particular heliocentric dis- 
tance is then easily calculated; Table I gives the mass-loss at r = 0.8 AU 
for various assumed r . 
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TABLE I. Bast Mass -Loss at 0.8 AU 


r: ' 

(AU) 

dM/dt 

/ 4 -1\ 

(10 g sec ) 

1.00 

5.93 

1.29 

5.47 

1.50 

5.17 

2.00 

4.78 

2.808 

4.42 


Thus, for this simple approximation, the mass-loss rate is not sensitive to 
the value of r , and dM/dt =- 5 x 10 4 g sec -1 is a typical value for r = 0.8 AU 

3. Velocity Distribution 

The evaporation of cometary ices (primarily H 2 0) produces a flow 
from the nucleus which accelerates dust particles to their terminal ejec- 
tion velocities within - 10 nuclear radii (Probstein, 1968; Delsemme and 

Miller, 1971) . Dels emm e and Miller (1971) calculate the largest particle 

, c -3 

which the expanding vapor can lift against gravity; for p - 0.4) g cm , 

R =1.7 km and p = lg cnf^ we have 
n n 

a =- 3.14 x 10 ^ Z cm 

rinoy 

-2 -1 

where the rate of evaporation, 2 (molecule cm sec ), depends on helio- 
centric distance. In Figure 3 we show the velocity distribution from Del- 
semme and Miller (1971) expressed in terms of a/a max . Using the evaporation 

17 

rate as reported by Marsden, Sekanina and Yeomans (197^)> at 0.8 AU 2 4.9x10 
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(KM/SEC) 



Figure 3 Terminal velocities of non-volatile particles* 





-2 

■aid a =“ 15.5 cm. Since Z varies roughly as r , larger particles can 
max 

\><i ejected at r <0.8 AU. Also, as can be seen iri Figure 3, at smaller 
heliocentric distances all particles tend to have nearly the same terminal 
velocity, while at larger distances v~ a (approximately). If, as is 

likely, Comet Enke is not covered with ices, then the effective Z(r) should 
be reduced, thereby also reducing a fflax and v(a). 

For calculating sunlight scattering from dust, we divide the particle 

sl/.e distribution into six ranges and assign appropriate mean radii to each 

O l/2 

range. The root-mean- square, (< a >) ' , in used to calculate the scatter- 

1 /? 2 

ing by each size range . While (< a ' >) is used for determining the terminal 

velocity, the result is not sensitive to the particular kind of averaging 

which is used because the velocity is a slowly- varying function of particle 

radius. The size ranges and derived properties are given in Table II, 

with m/dt = 5 x 10 4 g sec -1 and a = 15-5 cm as is appropriate for a hello- 
* max 

centric distance of 0.8 AU, 


TABLE II. Particle Sizes and Terminal Velocities (0.8 A U) 


Mange 

( | un) 

Production Rate 
(see -1 ) 

372 - 

(< a >) 

(tun) 

Geometric 
Scattering 
C ont rib ut ion( r /» ) 

v _ x 
(km sec ) 

l - 1.5 

1.51 

1 — l 
2 

* 

1.28 

9.8 

0.53 

1— 

1 

OJ 

1.79 

ln !4 
x 10 

2.04 

29.4 

0.49 

3 - 6.7 

3.68 

7 A 13 

x 10 

4.10 

24.5 

0.41 

6 . 7 - 14.4 

5.03 

12 

x 10 

9.35 

17.4 

0.33 

14.4-44.4 

8.77 

in 11 

x 10 

20.9 

15.2 

0.25 

*i4.4-a 

max 

2.34 

in 10 

x 10 

62.9 

3 . 65 

0.15 
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The first five size ranges used in the Mie scattering code are dis- 
cussed below; the final group contains particles sufficiently large that 
the geometrical approximation is adequate ("debris”). The scattering con- 
tribution (column 4 ) represents the relative total geometric cross section 
of each • size range. 

Because the particles attain their terminal velocities very near 
the' ‘comet ary nucleus, throughout the coma (where they will be observed) 
their spatial distribution is given by 

■ , , V, ' ; • • 

P. 3 

n (R) - ^ cm - 

1 . 4 jtR v. .• 

i 

where is the production rate and v_ L the terminal velocity for particles 

"t 7 Vl 

in the i ■ size range; R is the radial distance from the nucleus. 

This model shares with all current theoretical models the lack of 
asymmetry in ejection which must be present in order to produce the observed 
hon- gravitational forces. The size distribution of dust particles has the 
virtue of being related to observations of actual comets, although there has 
been no analysis of an old comet (such as Encke) using the procedures des- 

r 

cribed by Finson and Probstein (1968). Perhaps the lower curve in Figure 1 , 
giving Whipple T s ( 1968 ) distribution of interplanetary particles, would be 
a better representation for Comet Encke or the compilation by Dohnanyi 
(1972). The distribution used here contains a larger fraction of small 
particles (l < a < 5 pm) which are relatively efficient scatterers, so our 
results may represent an upper limit to the sunlight scattering for a given 
dM/dt. 
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IV. 


TRAJECTORY 


Comet Encke is a periodic comet with a period of 3.3 years. It shows 
greatly reduced activity after perihelion relative to before, and shows 
very little continuum radiation at any time, indicating a low dust content 
(Taylor et al. , 1973). Encke never shows a type II (dust) tail and in some 
apparitions has shown no tail at all. The observations from 1885-1991 indi- 
cates that Encke* s coma becomes observable at about 1.5 AU at the same time 
as the first appearance of its tail ( Vsekhsvyatskii, 19 64). The orbital 
characteristics of comet Encke are, given’ in Table 


TABLE III. Orbital Elements of Comet Encke 


Orbital period (year) 

3.30 

Aphelion distance (AU) 

4.09 

Perihelion distance (AU) 

0.34 

Orbital inclination (deg) 

12.4 

Velocity at 1 AU (km s 

37.1 

Velocity at perihelion (km s 

69.9 

Orbital eccentricity 

0.847 


5 

Vsekhsvyatskii (1964) has estimated a nominal value of 10 km for the dia- 
meter of Encke* s coma reduced to a Sun-comet distance of 1 AU. Observational 
data on tail length are given by Yoemans (1973) and TRW (1972). 

Many different kinds of missions to Encke have been proposed including 

fast and slow flybys. They do have certain elements in common. A nominal 

2 3 

miss distance is about 10 -10 km. One proposed probe trajectory calls for 

a rendevous at a heliocentric distance of about R = 0.8 AU. We have chosen 

s 

this distance to calculate the absolute,* and relative luminosities of the dust, 
debris, and nucleus. 
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V. 


MIE SCATTERING 


In this section the computer code used to calculate the Mie scatter- 
ing of sunlight from the cometary dust is described. Results are then presented 
for scattering of dust particles in the coma of Comet Encke. 

A. ' Mie Theory and Computer Code 

Particulate radiation effects in particulate clouds can be calculated 
provided that the optical constant of the radiating particles are known. 

Gal and Kirch ( 1973 and Gal, 1974) reported a new computer code called GMIE 

; • * 

to calculate scattering cross sections of interest. Calculations are based on 
MIE theory for a specified particulate column with a given complex index for re- 
fraction, range of wavelength, temperature, and particle sizes. 

It should be noted that the GMIE code can be used to obtain scattering 
cross sections or differential cross section for any given spherical particle or 
particulate column with a given size distribution and index of refraction. 

1. Electromagnetic Radiation Scattered by Absorbing Sphere 
The passage of electromagnetic radiation through a particle cloud is 
generally accompanied by removal of a fraction of the energy from the incident 
beam. This fraction may be partly absorbed within the particles and partly 
scattered -- i.e., reappear in the same direction as well as in other directions. 
The characteristics of the scattered radiation are determined by the wavelength 
X of the incident radiation, the refractive index (m = in^) of the particles 
and size as well as the shape of the discrete particles in the medium. For the 
radiation studies, we assume spherical particles with radius r. 
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The interaction of an electromagnetic wave with an absorbing sphere is described by 
the MIE theory (Fig. 4) and is discussed in detail elsewhere (van de Hulst, 1957) 
Classical MIE theory gives the intensity l(w/crn 2 ) at a distance H and angle B of the 
radiation scattered from a, single spherical particle of radius r exposed to parallel 
monochromatic radiation oi. intensity X^. 

J_ = « ( 0 ) ( 1 ) 

t 2(2 ir/A. f if 

where S ( 0 ) is the scattering function per unit particle. 

A computational scheme and computer code were provided by J. Dave of IBM (1968). 

These are based on the Ricatti-Bessel functions and give the MIE scattering cross 
section of absorbing spherical particles. 


Computations are valid for all values of size parameters that occur in the theory re 
gardlcss of whether they are large or small. The index of refraction of the particles 
are given in terms of their real and imaginary parts. 



Fig. 4 Definition of Single -Particle Scattering 

The extinction and scattering cross section arc calculated from various combinations 
of the sum and products of the coefficients and b R . The usual expressions for 

a and b are: 
n n 



r n ( y) - m yy> *;,<*> 
a » ; ^,’(y)T n (x) - nr^yTl^w 

m ! ' , n (y) ,f n l ' x) ‘ f„<y> (x) 

'’n " nTi^yTI^(x) - V~F)T|fW 


where. 

nv ~ the complex index of refraction 

n = a positive integer 

ip and 4 = the LUcatti-Bessel functions defined by 



ip (k) 

Y n 




J n+ 1/2 


(z) 


4 <*) 



u 


( 2 ) 

n+ 1/2 


<*) 


( 2 ) 


( 3 ) 


( 4 ) 

( 5 ) 


with J ^ 1/2 and H nU/2 the Bessel functions of lirst and third tyP es - The P rinie 
denotes differentiation with respect to the argument of the function and 

x = 2 7T r/A (6) 

y = mx (?) 


where r is the particle radius and A is the wavelength. 


Once these scattering coefficients are defined, then extinction and scatteiing nominal 
izod cross section can bo calculated with the following equations: 

oo 

Q,,XT ” i I (2 ” + JW < n n + V (8) 

x n - 1 
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q scatt 


uO 


-I I < 2 11 + 

X 1 

11= 1 




Scattering eiTioioncy (i.c., normalized scattering cross section) given by (.)) 
cmn also be written as a function ol scattering angle and the intensity in that direction: 


q scatt 


2 

x 


/ h <°> + *2 H 

o 


sin 0 dO 


( 10 ) 


who re 

i x (0) = |S T (0)| " 
i 2 («) - iS 2 (0)| 2 


and 


i 1 (0) + i 2 (0) = S(0) 


(11) 


i (0) and [< ( 0 ) refer, respectively, to the intensity of light polarized perpendicularly 
and parallel to the plane through the direction of propagation of the incident and scat- 
tered beam. These intensities are given in terms of the complex amplitude functions 
Sj (0) and S 2 <0) : 

a 7 T (cos 0) i- 1) r (cos 0)1 (l 2 ) 

n n n n i 


t and 

oO 

... V XilXJL 
s 2 (0) " Zj 11 ( n + 1) 

n = 1 



s (0) 


tL—j 


2 n \ 
n (n i 1) 
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The phase functions tt and appe 
terms of the Legendre polynomials, 


aring in Eqs. (12) and (12) can be expressed in 

P , as follows: 
n 


7 r (cos 0) 

n 


d P (cos 0) 
n 

d cos 0 


9 d7r T1 ( cos 

r n (COS 0) - COS 0 7T^ (COS 0) - si tP 0 d cos ~ 0 


(14) 


The detailed computational method to obtain these Legendre polynomials with recur- 
rence relationships is discussed elsewhere (Dave, 1968). 


If there is no absorption, be. , - 0 , then ^SCATT ‘ 

the efficiency factor for absorption, is given by 


Otherwise, » 


Q ABS q ext q scatt 


(15) 


Equation (0) or (10) gives the total scattered intensity. On the other hand, one may 
need to know the fraction of energy scattered into the forward and backward directions. 
This is obtained in terms of normalized cross section given in Eqs. (16) and (17), 
respectively 


and 


tt/2 


Q 


TUANS = \ I [*1 (0> + ‘2 (fl > 

X j 


sin 0 do 


(16) 


Q REF 


71 

= \ f f *1 <») + *2 < fi) 

x i/2 


sin 9 d0 


(17) 


2. Solar Scattering Cross Section for a Spherical Particle 

With the help of the basic MIE scattering theory, particular scattering due to 
sunshine will be treated through an equivalent normalized differential solar- 
scattering cross section. The sun is approximately half a degree as observed 
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from the earth, and calculations similar to those derived by GalandKirch (1973) 
for earthshine scattering would require evaluation of scattering function S(0) 
for smaller ’than 1 deg, which is our present minimum step size. Approximate 
differential cross section will be obtained by averaging the basic MIE scattering 
function over a small angular increment (A9) to account for the 1 irate solar disk 



Fig. 5 Solar Scattering Model 


From Fq. (1) 



(IB) 


The normalized dUL'erontiiU cross section can be obtained by dividing Eq. (18) with the 

2 

geometrical cross section G - r tt 


(dQ) 



(19) 


Equation (19) yields the 


normalized scattering cross 


section for a single spherical 


particle. 


3, Column Average Cross Section and Radiance 

a. Column Scattering 

The extension to the case of scattering from a column of particles is straight- 
forward upon adopting restrictive assumptions. First, the distance R must be 
large compared to the column dimensions so that scattering angle is essentially 
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constant for a given I q direction. Also, the column must, be* optically thin. 

This requirement means that: (l) each particle is exposed to the incident 
intensity I , regardless of location; and (2) the radiation scattered from 
one particle does not Interact with others as it passes from the column; i.e., 
there is no multiple scattering. 

These assumptions allow the scattered intensities to be added. With the help of Eq. (19) 
we can deline a cloud averaged differential cross section, o . 

a - r i 7r ( N / N t ) 2 » r " 1 ) ( 20 > 

i 

where 

N • = . the total number of particles 
t ' 

N. - the number i particles 
i li) 

r = the radius of the i particle 

i • 

b. Column Eadiance 

Figure 6 shows the geometry; the solid angle of the instrument is given by V = A/L , 
and the scattering cloud volume is V = AH, where A is the column surface area, and 
i is the depth of the column. The total number of particles in the scattering values is 



Figure 6 Definition of Scattering Values 
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N - nAI , whore » is the particle density. The incident intensity l^watt/em . *0 
can be calculated assuming the sun as a blackbody at a temperature I sun (>00 ° K * 


and sun is approximately half a degree as observed from the earth 


[ - Ik (T - 6000) 


/ -2 - 1 \ 
(A«) ( V/ cm pm ) 


(21) 


where 


B , 


1.19 


io 4 ( iA ) 5 


exp 


i. 4384 x 10 
X T 


4 


- 1 


. -2 -1 -Is 

( w cm sr pm ) 


( 22 ) 


whore A is tho wavelength In microns, and (An ) is solid angle of the sun viewed from 
the column (6.8 10 ^ s:r) at 1 AU. 


The column radiance is defined as 


% = -A 


, _2 -1 -Is 

( W cm sr pm ) 


(23) 


Substituting Eqs. (10 through (22) into Eq. (23), we obtain 


N, 


= B x n. r 1 (A«) S (dQ). rjir (N./N,.) (W cm' sr 


' 2 'tun' 1 ) (24) 


i = 1 


Equation (24) has been programmed into our GMIE code; additional inputs required are 
particle size, its distribution, and the column depth, i. 


26 



B. Mie Scattering From Cometary Past Particles 

We have inputted the comet Eneke dust distribution into the GMIE code 
for various indices of refraction. The normalized scattering and absorption 
cross sections given in Eqs. (lO)and (15) a re plotted in Figure 7 ^ or an index 
of refraction of I.7-.05 i. Figure 8 shows a plot of the scattering function 
as a function of scattering angle for a dielectric particle of O.65 pm radius, 
with an index of refraction n = 1.J-.05 i* More plots of particle scattering 
functions are found in Appendix A showing the effect of particle sizes and 
indices of refraction. Also, the angular dependence of the polarization is 
plotted for a few select cases. To calculate the scattered sunlight from the 
comet, we chose an index of refraction of I.7-.O5 i to represent cometary 
dust particles which are probably primarily silicates containing metallic 

elements , (Wickramasinghe and Krishna Swamy, 1968). Column, densities are calcu- 
lated in Appendix B. The results of the Mie scattering calculations are presented in 
Table IV where the wavelength dependence of the scattered radiance from Eq. (23) 
is presented as a function of the Sun -Comet -Probe angle. The functional dependence 
on wavelength and angle of the average differential cross section a from Eq. (20) 
is displayed graphically in Figs. 9 and 10 for a few selected cases. Figure 9 shows 
the effect of view angle on a at k = 0. 5 pm. The minimum scattering at this wave- 
length can be seen to occur in the 6-90° region. The wavelength dependence of 
the scattered radiance is displayed in Figure 10 for several view angles. It 
would appear that the scattering in the 0-90° range is relatively insensitive to 

wavelength dependence. 
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Figure 7 Scattering efficiency vs 0 Mie parameter for n=1.7-.05 i. 
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Figure 9 Angular dependence of average scattering cross section, 
A. = 0.5 dm 


3o 


cr (^m SR 



Figure 10 Wavelength dependence of average scattering cross section. 


31 




TABLE IV. Radiance of Comet Dust at R = 0.8 AU 

„ / -2 -1 -lv S 

N ( w cm sr pm ) 

K 


Sun-Comet -Probe Angle 



0° 

30° 

6o° 

90° 

120° 

150° 

H 

CC 

o 

o 

0.3 

2-7 -0 

2.4(-6) 

2.3( -6) 

2.2( -6) 

6.2(-6) 

1.2(-5) 

1.4 

0.4 

4.8(-6) 

4.2( -6) 

4.l(-6) 

4. l( -6) 

8.3(-6) 

2. 3( -5) 

1.5 

>v 0.5 

5.4(-6) • 

4.6(-6) 

4.4(-6) 

5.4(-6) 

8-9(-6) 

2.6(-5) 

1.1 

( fjjrri ) 0 . 6 

5.0( -6) 

4.l(-6) 

4.1( -6) 

6.o( - 6 ) 

9-0( -6) 

4.2(-5) 

6.8( -1) 

0.7 

4.3(-6) 

3.5(-6) 

3.3(-6) 

4.7(-6) 

6.8( -6) 

9.9(-6) 

4.2( -1) 

0.8 

3.7(-6) 

2.9(-6) 

2.5(-6) 

3.6( -6) 

6.3(-6) 

2.8(-5) 

2.7(-l) 

0.9 

2.9(-6) 

2.3( -6) 

2.2( -6) 

3.4( -6) 

4.0(-6) 

3.o(-5) 

1.8(-l) 

1.0 

2.6( -6) 

1.8(-6) 

1.8(-6) 

3.K-6) 

l-4(-6) 

l.8(-5) 

1.2(-1) 


Those results, which have been calculated for a heliocentric distance of 
R = 0.8 AU, can be approximately scaled to other distances by the law 

f7 * 

w. (0.8) , (25) 

s 

where Q accounts for the terminal velocity dependence on heliocentric dis- 
tance and R is in AU. Equation (25) reflects the inverse square dependencies 
s 

of the solar irradiance and dust production rates. The heliocentric dependence 
of the terminal velocities arises primarily in the case of Encke from the 
fact that the acceleration of the dust particles is dependent on the number 
of collisions the particles undergo with effusing gas molecules. Hence, 
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t'ho terminal velocities depend on particle sizes and the vaporization rate 
oi; gas from the nucleus. In the case of Encke dependence on the temperature 
of the nucleus is of lesser importance. For the particle sizes of importance 
here the heliocentric dependence can "be approximated by setting Q = 0. 1 !-. This 
law will hold for out to about 1.5 AU where the dust production ceases 
due to the threshold effect of ice evaporation. 
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VI. 


SCATTERING BY DEBRIS PARTICLES 


_2 

For particles in the 10 -10 cm range, scattering can be treated 

essentially as isotropic with a geometrical cross section. Thus, the 
differential cross section can be written as 


2 

, it a a 
da n 


dU 


4* 


(26) 


where the albedo a n is taken equal to that of the nucleus, 0.1. Calcu- 
lating an average differential cross section for R « 0.8 AU using the 

distribution in Table II results in 
__ 2 

a a 


da 

dQ 


n 


12-1 
9*9 * 10 pm sr 


(27) 


where a is the rms radius of large particles, a > 44 pm. The column den- 
sity for the debris particles is given by Eq. (B-3) which results in 



4jt 


1 

R . 
min 


-1 -2 

2,5 x 10 particles cm 


(28) 


Using Eqs. ( 27 ), (28) with Eq. (24), we then would get a scattered bright- 
ness for large particles at R r - 0.8 AU for X = 0.5 PM °f 

N = 8,9 x 10 ^ W cm ^sr \ari \ (29) 

K 


A.s can he seen from Table IV, this contribution is less than the small 
particle Mic scattering by a factor of 50» So it can be concluded, followin 
our model, that the small size dust particles contribute much more to the 
scattered brightness than the large size debris particles. 



VII. 


SCATTERING BY ICE PARTICIiES 


It. is difficult to precisely assess the ice particle distribution 
and production rates for Encke. We shall assume here that the ice particles 
can he represented by particles ■with a radius of 300 pm (Delsemme and Miller , 
1971). An upper limit for the ice production rate can he derived from the 
0G0-5 satellite observations of Lyman-a emission from Encke. The observed 
emission corresponded to a production rate of Q(H) - 5 x 10 2j atoms sr 1 

for hydrogen atoms at R = 0»715 AU (Bertaux et al. , 1973). Assuming that all the 

s 

H atoms came from ice particles results in a production rate for water Q(HgO) - 

27 ^ 

3 x 10 ' molecules s~ (Delsemme and Rud, 1973). For icy grains of 300 pm radius, 

1 9 -1 

we then have Q. ^ 2 x 10 s 
ice 

Using these numbers, it is possible to estimate the scattered sun- 
light assuming an albedo of 0.9. For the differential cross secti we 
have . 

- 5— -'2.0 x 10^ pm^ sr 1 , (30) 

and the column density would be 

~ 9 x 10”^ particles cm (3l) 

An estimate of the scattered light would then yield for X = 0.5 pm and 

R = 0.8 AU 

s 

N. = 6 x 10 ^ W cm ^ sr ^ pm \ (32) 

A, 

This result is of ihe same order of magnitude as the results for Mie 
scattering for cometary dust suggesting that scattering by ice particles 
surrounding the nucleus could be a major contributor to scattered sunlight. 

On the otherhand, our ice model is tentative and ice production rates could 
be. revised, downward. 


35 



VIII. 


SUNLIGHT REFLECTION FROM THE NUCLEUS 


We have assumed an albedo of 0.1 for the nucleus so the brightness 
of the nucleus following Eq. ( 24) is given by 

< uc = sr* < 33 > 

where AQ£ is the solid angle subtended by the sun at a given heliocentric 

-2 

distance, R - It varies as R . It is assumed here that the scattering 
s s 

is isotropic into 2 tc sr. For X - 0.5 iib and R = 0.8 AU, we have 

s 

Nf UC = 5=7 x 10" 3 W cm“ 2 sr" 1 um" 1 . (34) 

X 


IX. COMPARISON OF BRIGHTNESS OF NUCLEUS TO BACKGROUND BRIGHTNESS 

The contribution to the background brightness from the comet 
dust is tabulated in Table V. Taking one specific case we can compare 
the brightnesses from the nucleus to the background brightness. For the 
case with a heliocentric distance of 0.8 AU, X = 0.5 M-nL and view angle, 
a - 90°, the radiance levels are presented in Table IV. 


TABLE V. Radiances at R = 0.8 AU for X = 0.5 t-un, a - 90° 

s 


Nucleus 


, Particulate Scattering 
Dust Debris Ice 


Radiances 5*7 x 10 

(w c -i -is 

sr pm ) 


5=4 x 10~ 6 8.9 x 10“ 8 6 x 10“ 6 


The background from particulate scattering is 2 to 3 orders of magnitude 
below the intensity from the nucleus. The relative particle to nuclear 
brightness levels in Table IV are expected to vary less than an order of 
magnitude going into perhelion (0.34 AU) or out to "turn on" at about 2 AU. 
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Based on these calculations there should be no difficulty in viewing 


the nucleus from masking by particulate scattering. The particles are 
optically thin and the brightness levels are several orders of magnitude 
down from that of the nucleus. 



x. 


OPT IMIZATION OF IMAGING SYSTEM 


In terms of scientific value, resolution as high as possible is 
desired to examine the nature and composition of the surface. One of the 
problems with high resolution imagery and perhaps the limiting one is 
the large relative velocity of the spacecraft and comet which can be of the order of 
10 kin/sec. This, of course, would depend on the particular mission trajec- 
tory selected. Blurring of the image while a pixel cell is being exposed 
necessitates short exposure times. This requirement would argue for a 
high sensitivity-low noise system such as an image disector or a spin-scan 
camera. As an example of some numbers consider the case where the desired 
resolution element is ten meters squared on the surface of the nucleus and 

the bandwidth is 0.1 p. Then the photon rate corresponding to a resolution 

21 -1 -1 

element would be N - 1 x 10 photons- sec -sr . The solid angle subtended 


by a 5" diameter telescope at 1000 km would be about 2.5 x lo' 1 ^ ster 
and would result in photon rate at the detector of 3 x 10 6 photons/sec per 
pixel cell. If the transverse velocity of the comet with respect to the 


spacecraft is of the order of 10 km/sec then exposure times/pixel of no 
more than 1 msec are required. This corresponds to 6 x 10 J photons/pixel 
which is sufficient to detect intensity variations from the surface of 4-4 < 


It should be noted that the intrinsic 


readout noise of a vidicon system 


corresponds to about 1000 electrons /pixel. 


Other options which should be considered with an imaging system 
are use of polarizers to gain polarization information about the nuclear surface 
Because the illumination of the nucleus will be quasi-lambertian, viewing 
at diffeient phase angles will afford information regarding surface irregular- 


ities. Shadow heights can be used, to topographically map the surface. 
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of filters to block out line and band emission features in the bandwidth 
of the imaging system may be deemed wise. 

It is concluded that dust and debris scattering will not compro- 
mise the nucleus imaging mission and that high relative velocities warrant 
the use of a low noise imaging system such as a spin-scan camera. 
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APPENDIX A 


ANGULAR DISTRIBUTIONS FOR SINGLE 
PARTICLE MIE SCATTERING 
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APPENDIX B 


CALCULATION OF COLUMN DENSITIES 

The number of particles N passing through a spherical surface at a 
distance R from the center of the nucleus during a time At is given by 


N - S At = £ n^R) At 


4nR 


(B-l) 


■where P^, n^ and v^ are the production rate, the number density and the 
terminal velocity, respectively, of the particles in the ith size range. 


The number density as a function of R is then given by 


P. 

n (R) = — - - o 

1 v. 4*R 

x 


(B-2) 


The column density of particles n is obtained by integrating Eq. (B-2) 


from R . to R : 
mm max 


max P^dR 

Dco1 = i ^"r . 

man 1 


(B-3) 


R . is the radius of the nucleus (R ), and since v. is constant for 
mm n ' J i 

R > 10 R , the value of n .is insensitive to the value of R . Equation 
n- col max 

(B-3) is actually the lower limit to n since for R < 10 the particle 
velocities are less than The calculations of Delsemrae and Miller ( 1971 ), 

shown in their Figure 2, suggest that the particles are subject to 
approximately constant acceleration until they move at nearly terminal 
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velocity. The integral in equation (B-3) is most sensitive to the particle 
velocities at small R. As an upper limit to n col > suppose the actual 
particle velocities are given h y 


v(R) = 



R ^ R < 10 R 
n n 

R 10 R 

n 


(B-4) 


Here v q is the initial velocity at which particles are ejected from the 
nucleus. The column density for species i is then 


n 


col 


jits hr ^ + ' 0 . 9 ^ )]. 

n i 


(B-5) 


Taking ‘v = Oil v^ gives a column density 9«1 times the value obtained from 
equation (B-3) assuming v = v everywhere. We could, of course, go on to 
evaluate equation (B-3) using the assumption of constant acceleration to 
obtain somewhat better estimates for n ^ but the results depend on the 
unknown ejection velocity v . For the present, we simply point out this 
refinement which should be included in more comprehensive models. 


For calculating impacts suffered by a spacecraft passing near the 

comet, we have evaluated the column density for a linear path passing within 

R km of the nucleus (we will call this the impact parameter). The complica 
o 

t ions discussed above are not encountered because » R n * The spatial 
density is given by equation (B-2), and the column density is 


n col (a ’ R o ) = 


P(a) 


(B-6) 
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Since v(a) -* 0 as a -» a = the production rate P(a) must also go 
\ max 

to zero in order to avoid a singularity „ Since the terminal velocity for 
a > 10 a max ls given by 

v(a) “ 2.88 x 10 ' 3 (a/a^r^Cl-U/a^) 1 / 2 ] km/sec (B- 7 ) 

(from Delsemme and Miller, 1971), we cut off the production rate by a factoi 

[l-(a/a This does not affect the results presented in Section III, 

-3 

The production rate for a > 4.44 x 10 cm is given by 


P(a) = 7-298 x 10 M dugt a' 5 [l-(a/a n]ax ) 1/ ' 2 ] sec ' 1 (B- 8 ) 


° 4 -1 

where M, j =** 5 x 10 g sec and a = 15*5 cm at heliocentric distance 
dust max 

0,8 AU, while M. + 04 2,77 x Hr g sec and a = 98 » 2 cm at 0.34 AU 

ClUSo 

(perihelion) . 


Using equations (B-7) and B(- 8 ) in equation (B - 6 ), then integrating 

over particle sizes from a to a , we obtain the cumulative column density: 

max 

a 

max 

n(a,R 0 ) = J n co i da (B-9) 

a 

= 1.813 x 10' T [ a' 3 ' 5 - a^ 5 b (B-10) 

a R 
max o 

this represents the number of particles which have radii a or larger. 

Table B-I gives representative values. 
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TABLE B-I„ Column Densities For an 

Impact Parameter of 100 Ism 


Heliocentric 
Distance (AU) 


n(a, 100 kmj 
(cm" 2 ) 


-2 

a = 10 cm 

a = 10 ~ cm 

a = 1 cm 

0.34 

5.1 x 10' 3 

1.6 x 10 ^ 

5.1 x 10"^\ 

0.8 

2.3 x 10' 3 

7.3 x 10“ T 

2.3 x 10- 10 


For given values of a and R q , n varies approximately inverse^ as the 
heliocentric distance because M dust and a fflay both vary approximately as 
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PARTICLE RADIUS (CM) 


Figure B1 Cumulative column density vs. particle radius. 
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